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Identification and expression analysis of AP2/ERF transcription 
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ERF 分 为 ERF 和 DREB 两 个 亚 类 ，ERF 包含 ERF 
5 类 ， 这 与 模式 植物 拟 南 芥 AP2/ERF 的 
橄榄 的 AP2/ERF 家 族 
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， 这 与 转录 组 分 析 结 果 一 致 。 
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| Abstract: In order to explore the response mechanism of AP2/ERF gene family in the water stress of Olea 


europaea, this study performed transcriptome sequencing on the roots and leaves of two cultivars 'Frantoio' and 


'TYZ-1' that were under drought and flooding stress. And based on the whole genome data, the protein 
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physicochemical properties, gene structure and system evolution of AP2/ERF transcription factor in O. europaea 
were analyzed. At the same time, the difference in gene expression of AP2/ERF transcription factor related to 
water stress in the two O. europaea cultivars was analyzed by transcriptome sequencing data and verified by 
fluorescence quantitative PCR. The results were as follows: (1) A total of 110 AP2/ERF gene family members 
were identified in O. europaea. The amino acid size of the 110 proteins is 173-717bp, there is no signal peptide 
and it is a non-secreted protein. The phylogenetic tree was constructed between O. europaea AP2/ERF and model 
plant Arabidopsis AP2/ERF protein. It was found that O. europaea AP2/ERF protein was divided into four 
categories: AP2, RAV, ERF and Solosist. Among them, ERF is divided into two subtypes: ERF and DREB. ERF 
includes six subtypes of ERF B1 to ERF B6, and DREB includes six subtypes of DREB Al to DREB A6, which 
is consistent with the classification of the model plant Arabidopsis AP2/ERF. Each subfamily contains AP2/ERF 
proteins of O. europaea and Arabidopsis at the same time, indicating that the AP2/ERF family of Arabidopsis and 
O. europaea are similar in evolution. (2) The analysis of gene structure and conserved elements found that the 
proteins of the same subfamily of O. europaea AP2/ERF have the same gene structure and conserved elements. 
Combining gene expression with genes with known water regulation functions in the evolutionary tree, it is 
preliminarily speculated that OeAP2-75, OeAP2-97, OeAP2-101, OeAP2-23, OeAP2-13 are closely related to the 
water regulation of O. europaea, OeAP2-13, OeAP2-28, OeAP2-104, OeAP2-75, OeAP2-80, OeAP2-50 have 
different expression levels in the two varieties. It is speculated that this may be the reason for the different 
resistance of  'Frantoio' and 'TYZ-1'. The RT-qPCR technique was used to detect the expression changes of O. 
europaea AP2/ERF gene under different stresses. The results showed that OeAP2-101, OeAP2-28, and OeAP2-42 
were significantly up-regulated by water stress, which was consistent with the results of transcriptome analysis. 
The results of this study can lay a foundation for the research on the stress resistance expression and gene 
function of the AP2/ERF family genes of O. europaea, and provide a method and theoretical basis for the 
selection of drought-resistant and flood-tolerant rootstock varieties of O. europaea. 

Key words: Olea europaea, AP2/ERF transcription factor, bioinformatics, gene expression profile, water stress 


植物 在 自然 界 生 长 发 育 和 进化 的 过 程 中 ， 通 常 面临 着 诸如 干旱 、 水 济 、 极 端 高 低温 、 高 盐 等 的 非 生物 
胁迫 ， 这 对 植物 的 生长 、 发 育 产生 了 负面 影响 。 因 此 植物 通过 多 种 信号 转 导 途径 使 转录 因子 与 下 游 基因 的 
顺 式 作 用 元 件 相 结合 ， 调 控 下 游 防卫 反应 基因 的 表达 ， 从 而 在 细胞 和 分 子 水 平 上 对 环境 胁迫 作出 反应 ( 张 
XC, 2016)。 目 前 研究 发 现在 非 生物 胁迫 中 起 重要 作用 的 转录 因子 包括 AP2/ERF、WRKY、NAC、MYB、 
ZFP, bHLH 等 家 族 (Xu et al., 2011). 

AP2/ERF 是 植物 中 最 大 的 转录 因子 家 族 之 一 ， 涉 及 植物 的 生长 发 育 及 各 个 生理 过 程 ， 参 与 植物 非 生 
物 胁迫 调控 机 制 ,AP2/ERF 家 族 蛋 白 的 主要 特征 是 含有 1 个 或 2 个 AP2 结合 域 ,每 个 AP2 结合 域 约 含 58~70 
个 氨基 酸 残 基 (Cao et al., 2020)。 根 据 AP2 家 族 和 蛋白 结 构 域 中 特征 元 件 的 种 类 和 数量 不 同 可 将 其 分 为 AP2、 
ERF, RAV, Soloist 4 个 亚 家 族 (Wu et al., 2015)， 各 个 亚 家 族 在 植物 体内 有 不 同 的 功能 ，AP2 亚 家 族 包含 
2 个 AP2 保守 结构 域 ， 主 要 与 植物 生长 发 育 及 细胞 生长 分 化 有 关 ( 纪 晴 等 ，2018)，RAYV 含有 1 个 AP2 Zi 
构 域 和 1 B3 DNA 结合 域 ,通常 在 乙烯 、 蔚 蔡 素 内 酯 及 一 些 生物 和 非 生 物 胁 迫 中 起 作用 ( 柯 希 望 等 , 2020); 
ERF 和 Soloist 只 包含 1 个 AP2 结构 域 ， 其 中 ，ERF 又 分 为 DREB 和 ERF 两 个 亚 类 ， 通 常 与 植物 的 生物 
胁迫 、 干 旱 、 高 盐 、 低 温 、 热 胁迫 、 多 重 胁迫 相关 ( 刘 志 薇 等 ，2014; 苟 艳丽 等 ，2020)。 大 量 研究 证 明了 
AP2/ERF 转录 因子 在 植物 生长 发 育 及 非 生 物 胁迫 中 起 关键 的 调控 作用 ， 例 如 沙 柳 (Salix cheilophila) 
SpsDREBS 基因 在 干旱 胁迫 下 表达 下 调 ( 王 雷 等 ，2021)， 拟 南 芥 Cdrapidopsis thaliana) F EA H & (Ipomoea 
batatas) 的 IbRAP2-12 基因 提高 其 植株 的 耐 盐 性 和 抗旱 性 (Li et al., 2019)， 过 表达 OsERF71 Hei KAR (Oryza 


LEH 


Ne 


pap 


sativa) 的 耐 旱 性 (Ahn et al., 2017). 。 随 着 近年 植物 基因 组 数据 相继 公布 ， 发 现 不 同 植物 中 AP2/ERF 转录 因 
子 的 数量 不 同 ， 比 如 中 国 樱花 (Cerasvs serrulata) 中 鉴定 出 了 68 个 AP2/ERF 转录 因子 家 族 成 员 (Zhu et al., 
2021), "Sz*(Dactylis glomerata) d FATE(Betula platyphylla Suh) 中 分 别 有 193 个 和 45 个 AP2/ERF 转录 因子 
家 族 成 员 CXu et al, 2020; 张 文 慧 等 ，2020)， 但 目前 对 油橄榄 (Oiea exropaea) 基 因 组 中 的 AP2/ERF 转录 因 
子 的 研究 未 见报 道 。 
橄榄 为 木犀 科 (exropaew) 木 犀 榄 属 (Oieg) 油 料 作 物 ， 原 产 于 地 中 海 沿 岸 ， 适 应 夏季 长 、 热 、 干 燥 ， 冬 
季 温 和 多 雨 的 半 干 旱 气 候 (Amira et al., 2020), 是 世界 著名 的 亚热带 果树 和 重要 经 济 林木 ( 牛 二 利 等 , 2021)， 
其 全 身 是 宝 ， 橄 榄 叶 和 橄榄 油 中 富 含 多 种 有 益 活性 物质 ， 易 于 人 体 吸 收 ( 程 子 彰 等 ，2014)， 橄 榄 油 在 西方 
被 誉 为 “液体 黄金 ”, 是 世界 上 唯一 直接 采用 鲜果 冷 榨 工 艺 以 自然 形态 榨取 的 木 本 植物 油 ( 邓 从 静 等 , 2011). 
随 着 社会 经 济 发 展 , 橄榄 油 消费 需求 急剧 增长 , 橄榄 油 市 场 长 期 供不应求 、 价 格 居 高 不 下 ( 赵 梦 炯 等 , 2021). 
云南 油橄榄 栽培 由 于 引种 地 和 原 产 地 生态 条 件 差异 还 面临 很 多 问题 , 夏季 湿润 多 雨 , 冬季 干旱 的 气候 与 原 
产地 气候 正好 相反 ， 因 此 选 育 夏季 耐 洲 、 冬 季 抗 旱 的 油橄榄 砧木 品种 是 解决 云南 油橄榄 裁 培 问题 的 关键 。 
该 研究 基于 课题 组 前 期 田间 试验 结果 田园 1 号 (TYZ-1 号 ) 在 水 胁迫 下 抗 道 性 强 于 ' 佛 奥 ， 为 探究 2 个 品种 
抗 逆 性 差异 的 原因 ， 该 文 对 油橄榄 AP2/ERF 基因 进行 挖掘 鉴定 ， 并 分 析 其 蛋白 理化 性 质 、 系 统 发 育 及 水 
胁迫 下 的 基因 表达 模式 ， 为 进一步 研究 油橄榄 AP2/ERF 家 族 基 因 的 抗 着 性 表达 及 功能 基因 的 挖掘 黄 定 基 
as 础 ， 也 为 油橄榄 选 育 抗旱 、 耐 涝 砧木 品种 提供 了 方法 和 理论 依据 。 

O 研究 发 现 油橄榄 AP2/ERF 基因 家 族 包含 110 个 成 员 ， 分 类 与 模式 植物 拟 南 芥 一 致 ， 其 在 进化 水 平 上 
定 的 相似 性 。AP2/ERF 家 族 中 转录 因子 的 数量 不 但 与 物种 的 基因 组 有 关 ， 还 与 植物 长 期 进化 过 程 中 
CN ”外 界 环境 的 影响 有 关 ， 与 枣 (Ziziphus jujuba ) 基 因 家 族 中 145 个 ( 李 继 东 等 ，2020)， 银 杏 (Ginkgo biloba) + 
Liu 61 ARAR, 2020), (Salix arbutifolia) P 173 ^^ (Rao et al., 2015)AP2/ERF 家 族 成 员 相 比 ， 油 橄榄 
中 AP2/ERF 转录 因子 的 数量 属于 中 等 类 型 ， 这 一 结果 推测 油橄榄 在 自然 进化 过 程 中 ，AP2/ERF 基因 家 族 
受到 一 定 的 外 部 环境 压力 。 


1 材料 与 方法 


1.1 试验 材料 获取 及 转录 组 测序 

c 试验 材料 为 油橄榄 'TYZ-1 号 ?和 ' 佛 奥 : 一 年 生 插 揪 苗 ， 由 云南 省 林业 和 草原 科学 院 树 木 园 提供 ('TYZ-1 
号 是 绿化 用 的 油橄榄 种 子 通过 杂 合 培育 出 的 实生 苗 中 选 育 的 优良 无 性 系 ， 其 耐 兰 薄 、 抗 旱 耐 涝 性 强 。 ' 佛 
奥 ? 是 我 国 通 过 国家 审定 的 油 橄 槛 恨 种 ， 中 国 油 橄榄 适 生 区 大 都 有 种 植 ， 但 抗 逆 性 不 强 )。 试 验 前 期 选取 2 
个 品种 油橄榄 各 30 株 幼 苗 分 干旱 、 水 淹 及 对 照 组 于 2019 年 9 月 3 日 开始 , 对 2 个 品种 插 捅 苗 进 行 干旱 、 
水 淹 及 正常 培养 处 理 ，2019 ik A 月 8 日 采集 不 同 处 理 (FPE KE RE 下 油橄榄 (‘TYZ-1 5 R) 
TH Br EORR 3 个 混合 样 。 所 采集 的 样品 当即 处 理 后 用 液 氮 速 冻 ，-80'C 下 保存 ， 按 照 RNA 试剂 盒 
(Qiagem 说 明 提 取 RNA, bh RNA A 为 模板 合成 cDNA， 并 构建 文库 。 委 托 上 海派 森 诺 生物 科技 有 限 公 司 利 
用 二 代 测 序 技术 (Next-Generation Sequencing, NGS) Æ F Illumina 测序 平台 对 样品 文库 进行 双 末 端 
(Paired-end，2*150bp，PE) 测 序 。 
1.2 OeAP2/ERF 蛋白 序列 的 获取 
油橄榄 的 全 基因 组 数据 从 NCBI 上 下 载 ， 搜 索 获 得 油橄榄 AP2/ERF 转录 因子 的 cDNA. 序列 及 和 蛋白 序 
列 ， 利 用 NCBI 进行 BLAST 同 源 序列 比 对 及 和 蛋白 结构 域 分 机， 除去 无 AP2/ERF 结构 域 的 蛋白 序列 。 
1.3 OeAP2/ERF 蛋白 理化 性 质 分 析 

利用 ExPASy protparam tool 在 线 软件 分 析 AP2/ERF 蛋白 理化 性 质 ， 通 过 signalP-5.0 和 Euk-mPLoc 
2.0 server 预测 蛋白 的 信号 肽 和 亚 细 胞 定位 ， 利 用 Prabi 在 线 软件 (NPS@ : SOPMA secondary structure 


ne 


prediction) 预 测 重 白 的 二 级 结构 。 


1.4 OeAP2/ERF 蛋白 进化 树 构 建 及 保守 元 件 和 基因 结构 分 析 


利用 MEGA X 以 及 Figtree KK A FY 
行 参数 为 None。 油 橄榄 AP2/ERF KIKE A 


分 类 参考 拟 南 芥 AP2/ERF 家 族 蛋 白 的 分 


橄榄 和 拟 南 芥 AP2/ERF 蛋白 序列 构建 进化 树 构建 及 修饰 ， 执 


类 方法 ， 拟 南 芥 蛋白 


序列 从 拟 南 芥 基 因数 据 库 下 载 (https:/www.arabidopsis.org/download/index-auto.jsp? dir=%2Fdownload_ 


files% 2FProteins)， 利 用 MEME (http://meme.nbcr.net/meme/tools/ meme) 对 110 个 OeAP2/ERF 4z 
F: 基 序 位 点 分 布 情况 ， 选 择 重 复 次 数 不 限 制 ; 


构 域 进行 预测 ， 上 有 具体 参数 设置 如 
择 6， 其 他 参数 均 采 用 默认 值 。 


1.5 OeAP2/ERF 水 分 胁迫 基因 表达 谱 分 析 


通过 有 参 转录 组 测序 获得 ; 


FPKM 值 数据 ， 对 其 进行 排序 ， 
变化 最 大 的 15 个 基因 。 
1.6 实时 定量 荧光 PCR 


表达 交互 热 图 。 非 生物 胁迫 基因 表达 谱 


数据 的 处 理 及 筛选 采用 FPKM 对 表达 量 进行 标准 
删除 FPKM 无 变化 或 者 变化 范围 不 超过 


| 橄榄 AP2/ERF 转录 因子 基因 表达 谱 ， 利 用 派 森 诺 基 因 


保守 


为 了 验证 转录 组 测序 的 结果 及 目 
| 橄榄 OeAP2 基因 在 


荧光 定量 PCR 检测 ; 


(invitrogem) 检 测 特异 引物 的 PCR 产 4 
上 下 游 引物 (10 um/1)0.5 pL， 模板 (cDNA)1 uL, ddH20 10.5 uL. fi 
系统 , Foster City, CA, USA). PCR 反应 程序 : 变性 程序 (95 C, 10 min), WKE 
| 线程 序 (60 CE 95 C, 
正常 表达 的 内 部 调控 因子 ， 
生物 学 重复 和 每 个 生物 学 重复 的 3 个 技术 重复 。 
橄榄 荧光 定量 检测 引物 序列 


Table 1 Sequence of the primers for fluorescent quantitative detection of Olea europaea 


57°C, 10s; 72 C,15 s; 单 次 荧光 
后 冷却 至 40°C. V tubulin 基因 


标 基因 


的 表达 模式 ， 设 计 特 异 引 物 ( 表 
“TYZ-1 5 RI fh BL? f [ei] Ab H 
7]. 25 ul 反应 体系 的 可 选 参数 如 


5 倍 的 数据 ， 


RH 


测量 )， 熔 化 


作为 基因 


xli 


保守 结 


FE 基 序 的 数目 限制 选 


云 在 线 软件 作 基因 
化 ， 保 留 所 有 
筛选 出 基因 表达 量 


1), tubulin 作为 内 参 基因 ， 经 
下 的 表达 情况 。 用 SYBR Green 

F: 2x SYBR 绿色 主 混 
] PCR 热 循环 仪 (ABIi 7300; 应 


合 物 12.5 uL, 


JEW 


EH 45 (95C, 15 s; 
Da dl C, ERIM), HK 
通过 qPCR 分 析 每 个 样品 3 个 独立 的 


引物 名 称 Primer name 


序列 sequence 5'-3' 


OeAP2-101F CAACAATTCCAGCTCCTCCA 
OeAP2-101R CTGCTTCTTCAGCTGTGTT 
OeAP2-55F TACAGAGGAATCCGCATGAG 
OeAP2-55R CGACAATTCATGACATACTT 
OeAP2-42F TCAAGGTGGATCTGTGAGGA 
OeAP2-42R ATCACTGCTTATATATCGCA 
OeAP2-28F TCAGATATATGCCAACCATC 
OeAP2-28R TGACCTACTTGGGCTCGTGC 
OeAP2-23F GAGGAATATGTACAGAGGCA 
OeAP2-23R AGTCGAAGGTGGCGGATTGT 
tubulinF AGATGATGATGACGAGTATG 
tubulinR GACGATGTATCCGTCTATCA 
结果 与 分 析 


2.1 有 参 转录 组 测序 数据 整理 


油橄榄 'TYZ-1 号 及 ' 佛 奥 : 在 和 干旱、 水潭 、 正 常 浇 水 下 的 根 与 叶 共 12 组 文库 经 上 机 测序 ， 得 到 图 像 文 
件 ， 由 了 umina 测序 平台 自 带 软件 进行 转化 ， 生 成 FASTQ 的 原始 数据 (Raw Data)， 对 每 个 样品 的 下 机 数 
据 (Raw Data) 分 别 进 行 统计 ， 包 括 样 品名 、Q30、 模 糊 碱 基 所 占 百分比 、 以 及 Q20(%) 和 Q30(%)， 统 计 结 
果 见 表 1。 测 序数 据 结果 显示 ， 所 有 样品 Reads 总 数 在 38 191 054~46 939 278 之 间 ; 碱 基 总 数 在 728658 
100 bp 以 上 ;Q20 碱 基 百 分 比 均 大 于 96.91 96; Q30 碱 基 百 分 比 均 在 92.2 % 以 上 ; 模糊 碱 基 所 占 比 例 在 0.000 
769 % 以 下 ， 后 经 数据 过 滤 去 除 这 些 带 接头 、 低 质量 的 Reads. 

表 2 测序 数据 统计 


Table 2 Sequencing data statistics 


样品 名 Reads 总 数 ERB — 碱 基准 确 率 99% 模糊 碱 基 碱 基准 确 率 ” 碱 基 准确 率 
Sample Total Bases (bp) ”以 上 碱 基 总 数 ”百分比 99 % 以 上 的 碱 99.9 % 以 上 的 
number of Q30 (bp) N (%) 基 所 占 百分比 碱 基 所 占 百 分 
Reads Q20 (%) 比 Q30 (%) 
FDL 45048662 6791084042 6304559475 0.000268 97.36 92.83 
FDR 46939278 7074531496 6522903901 0.000291 96.91 92.2 
ODL 44538262 6725277562 6240276298 0.00025 97.29 92.78 
ODR 45559988 6866263812 6337361057 0.000 283 96.94 92.29 
FSR 38716652 5 807 497 800 5481226288 0.000556 97.97 94.38 
FSL 42 434914 6365237100 6014466495 0.000 75 98.02 94.48 
OSR 43986834 6598025100 6247047129 0.000768 98.08 94.68 
OSL 41168206 6175230900 5834000945 0.000749 98.03 94.47 
FR 41814712 6272206800 5901059489 0.000758 97.83 94.08 
FL 42 041 262 6306189300 5982899942 0.000764 98.21 94.87 
OR 38191054 5728658100 5429155675 0.00 076 98.11 94.77 
OL 41 495 232 6224284800 5912636170 0.000769 98.24 94.99 


F 表示 ' 佛 奥 '，O 表示 “TYZ-1 5’; D 表示 干旱 ; S 表示 水 淹 ; R 表示 根 ; L 表示 叶片 。 
F stands for 'Frantoio'; O stands for 'TYZ-1'; D is for drought; S is for flooding ;R is the root; L is the leaf. 


2.2 样品 基因 差异 表达 分 析 

采用 DESeq 对 基因 表达 进行 差异 分 析 ， 表 达 差 异 倍数 |log2FoldChange| > 1， 显 著 性 p< 0.05。 不 同 条 
件 下 的 样品 进行 两 两 比较 , 图 1 为 干旱 胁迫 下 油橄榄 2 个 品种 基因 表达 差异 分 析 结 果 簇 状 柱 形 图 ， 图 2 为 
水 淹 胁 迫 下 油橄榄 2 个 品种 基因 表达 差异 分 析 结 果 簇 状 柱 形 图 , 图 1 中 FDR 与 FDL 之 间 差 异 基因 数量 最 
少 ， 差 异 基因 总 数 为 92 个 ， 上 调 差 异 基 因数 为 45 个 ， 下 调 差异 基因 数 为 47 个 ; FR 与 OL 之 间 差 异 基因 
数 最 大 ， 这 是 两 个 品种 在 正常 浇 水 下 的 根 与 叶 之 间 的 基因 表达 差异 分 析 ， 差 异 基因 总 数 为 2 803 个 ， 上 调 
差异 基因 数 有 1 076 个 , 下 调 差 异 基因 数 有 1 727 +. 图 2 F OSR 与 OR 之 间 差 异 基因 数 最 少 为 1 550 个 ， 
上 调 差 异 基因 数 有 693 个 ， 下 调 差 异 基因 数 有 857 个 ; 除 FR 与 OL 外 ，FR 与 FL 之 间 差 异 基 因数 最 大 为 
2 791 个 ， 上 调 差 异 基因 数 为 976 个 ， 下 调 差异 基因 数 为 1 815 个 。 由 上 述 数 据 得 知 ， 在 干旱 和 水 淹 胁 迫 
下 油橄榄 不 同 品种 中 根 与 叶 之 间 的 差异 基因 数量 较 明 显 。 
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F RAN HER’; O 表示 “TYZ-1 Su D 表示 干旱 ，R 表示 根 ; 工 表示 叶片 。 下 同 。 
F stands for 'Frantoio'; O stands for 'TYZ-1'; D is for drought; R is the root; L is the leaf. The same below. 
1.ODL vs FR; 2. FDL vs ODR; 3. OL vs FDL; 4. FDL vs OR; 5. FDL vs ODL; 6. FDR vs FR; 7. OR vs OL; 8. 
ODR vs FL; 9. FDR vs ODR; 10. OL vs ODL; 11. FR vs FL; 12. FDR vs FL; 13. ODR vs FR; 14. ODR vs ODL; 
15. FL vs OR; 16. FL vs OL; 17. FR vs OR; 18. FDR vs ODL; 19. FDL vs OL; 20. ODR vs OL; 21. FDL vs FR; 
22. OR vs FDR; 23. ODL vs OR; 24. ODL vs FL; 25. FDR vs OL; 26. FDL vs FL; 27. ODR vs OR; 28. OR vs 
ODR; 29. FR vs OL; 30. FDR vs OR; 31. ODL vs OL; 32. FDR vs FDL 

1 干旱 胁迫 下 基因 表达 差异 分 析 结果 统计 


Fig. 1 Statistical diagram of drought stress gene expression difference analysis results 
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S 表示 水 淹 S is for flooding 


1. OSR vs OSL; 2. FSR vs OSL; 3. FSL vs OSL; 4. FSL vs OSR; 5. FSL vs FR; 6. FSR vs OSR; 7. OR vs OL; 8. 
OSL vs OR; 9. FSL vs OL; 10. FSR vs OR; 11. FR vs FL; 12. OSL vs OL; 13. FSR vs FSL; 14. OSR vs FL; 15. 


FR vs OL; 16. FSR vs FL; 17. FSR vs OL; 18. OSR vs OL; 19. FSR vs FR; 20. FL vs OR; 21. FSL vs FL; 22. 


FSL vs OR; 23. OSR vs FR; 24. OSL vs FL; 25. OSR vs OR; 26. FL vs OL; 27. FR vs OR; 28. OSL vs FR 


2 水 淹 胁 迫 下 基 


因 表 达 差 异 分 析 结 果 统 计 图 


Fig. 2 Statistical diagram of water-out stress gene expression difference analysis results 


2.3 油橄榄 AP2/ERF 转录 因子 家 族 成 员 鉴 定 及 理化 性 质 分 析 
究 从 油橄榄 全 基因 组 中 搜索 AP2/ERF 基因 ， 共 鉴定 出 110 个 成 员 ， 该 110 个 蛋 


AST 


白质 所 含 氨基 酸 


大 小 为 173-717bp， 平 均 氨 基 酸 个 数 为 344 个 。 通 过 ExPASy protparam tool 分 析 OeAP2/ERF 蛋白 理化 性 


Ji, 110 个 OQeAP2/ERF 和 蛋白 中 37 个 蛋白 的 PI 大 于 7 4 
OeAP2-6、OeAP2-28、OeAP2-109) 不 稳定 系数 小 于 40， 为 稳定 蛋白 


系 >40); 110 个 OeAP2/ERF 蛋白 的 脂肪 族 系数 均 小 于 100; Fil) ae 


白 都 不 存在 信号 肽 ， 为 非 分 泌 和 蛋白 ; 


5 细胞 定位 中 有 56 个 OeAP2/ERF 4 


于 细胞 质 中 ，25 个 同时 定位 于 旨 


FH 


胞 质 和 细胞 核 中 ，OeAP2-85 定位 于 细胞 质 和 线粒体 


EF 均 PI 为 6.63; 有 4 个 AP2/ERF 转录 因 
， 其 余 106 个 为 不 稳定 蛋白 (不 稳定 
Aas hk, RIL 110 个 OeAP2/ERF 


"f (OeAP2-97. 


AE 


白 定位 于 细胞 核 


现 OeAP2/ERF 蛋白 无 规则 卷曲 为 主要 结构 ，o- 螺 旋 为 次 要 结构 。 


2.4 OeAP2/ERF 蛋白 进化 分 析 与 分 类 


从 拟 南 芥 基 因 库 下 载 拟 南 芥 AP2/ERF 转录 因 
类 ， 从 每 个 大 类 中 挑选 1-12 个 AP2/ERF 转录 因子 作为 参照 。 利 用 MEGA XX 对 
白 及 拟 南 芥 32 AP2/ERF 4& 


KATE 
体 分 为 四 类 ， 其 
另外 两 个 成 员 (OeAP2-28、OeAP2-87) 
xE, m 
独 的 分 支 OeAP2-60, 


和 OeAP2/ERF 蛋白 ， 说 明 拟 南 芥 和 油橄榄 的 AP2/ERF 家 族 在 进化 水 平 上 有 一 定 的 相似 性 


DREB A4 


于 RAV 类 ; 拟 南 芥 中 Solosist 类 和 蛋白 与 RAV 
推测 其 属于 油橄榄 
子 家 族 其 他 成 员 差 距 较 大 ; 其 余 的 86 个 成 员 都 包含 1 AP2 

包含 ERF 和 DREB 两 个 亚 类 , 分 别 与 拟 南 芥 ERF 和 DREB RA 
B1~ERF B6, DREB 分 为 DREB A1~DREB A6, 除了 Solosist 外 , 其 余 每 个 亚 组 中 都 同时 包 


包含 1 个 AP2 £i sug 


子 ， 根 据 Wu 等 (201$) 对 拟 南 芥 AP2/ERF 转录 因 


白 进行 系统 进化 树 构 建 及 分 析 , 结果 如 图 3 所 示 ，AP2/ERF f& 
中 与 拟 南 芥 AP2 类 聚 在 一 起 包含 两 个 AP2 保守 结构 域 的 成 员 共 有 21 个 ， 被 分 为 AP2 类 ; 
HER, SMA RAV 聚 在 一 个 分 
Il Solosist 385974 — ^ 
sist Ej AP2/ERF 转录 因 


I 1^ B3 
KREE, ME RAV RA 
Solosist 类 成 员 ， 单 独 分 支 可 能 是 Solo 


AE 


呆 守 结构 域 ， 同 属于 ERF 28, ERF 类 


o 


AE 
E: 


，28 个 定位 
1; 预测 二 级 结构 发 


子 的 分 
橄榄 110 个 AP2/ERF 转 
白 总 


1M 


E 一 起 , 总 共 分 为 12 个 亚 组 , ERF 分 为 ERF 
& T AtAP2/ERF 


图 3 油橄榄 与 拟 南 芥 AP2/ERF 家 族 系统 进化 树 


ChinaXiv 


Fig. 3 Phylogenetic tree of the AP2/ERF family of Olea europaea and Arabidopsis 


2.5 油橄榄 AP2/ERF 保守 元 件 及 基因 结构 分 析 


利用 MEME 对 OeAP2/ERF 进行 保守 元 件 分 析 , 得 到 相关 性 最 高 的 6 个 保守 元 件 如 图 


OeAP2/ERF 都 包含 motif 1， 说 明 motif 1 是 油橄榄 AP2/ERF 的 保守 结构 域 ，y 


橄榄 AP2/E 


4 Bias, 所 有 的 
RF 转录 因子 的 


保守 元 件 及 基因 结构 还 与 其 分 类 相关 ， 在 油橄榄 AP2/ERF 分 类 中 的 AP2 类 转录 因子 100% 包 含 motif 1, 


96.4% 包 含 motif 2，89.7% 包 含 motif 3，55.2% 包 含 motif 4，37.9% 包 含 motif 5; 


可 能 是 由 于 在 进化 过 程 中 发 生 了 变异 ,导致 该 类 转录 因子 与 其 他 成 员 之 间 没 有 


包含 motif 1, motif 2 和 motif 3; DREB 类 均 包含 motif 1, motif 2 和 motif 6， 每 一 亚 类 中 的 保守 基 序 位 置 


RAV 类 只 


包含 motif 1, 


ET AY) EST d 


\ 


大 臻 相同， 基因 结构 也 一 致 ， 说 明 油 橄榄 OeAP2/ERF 在 进化 过 程 中 的 高 度 保守 性 


EFE. ERF 类 均 
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图 4 油橄榄 AP2/ERF 基因 家 族 保守 元 件 和 基因 结构 分 析 

Fig.4 Analysis of conserved elements and gene structure of Olea europaea AP2/ERF gene family 
2.6 油橄榄 AP2/ERF 水 分 胁迫 基因 表达 谱 分 析 
通过 基因 表达 谱 作 出 水 胁迫 基因 表达 交互 热 图 (图 5), 油橄榄 AP2/ERF 家 族 成 员 中 共 检 测 到 100 个 基 
因 表达 ,未 检测 到 Oe4P2-29、OeAP2-30、Oe4P2-40、OeAP2-56、OeAP2-59、Oe4AP2-72、OeAP2-81、OeAP2-87、 
OeAP2-88. OeAP2-105 基因 的 表达 ,推测 这 10 个 基因 未 参与 水 胁迫 响应 。 干旱 胁迫 下 , 与 正常 浇 水 相 比 ， 
“ 佛 奥 叶片 中 40% 基 因 表 达 上 调 ，27% 的 基因 无 表达 ，33% 的 基因 表达 下 调 ， 根 部 33% 基 因 表 达 上 调 ，9% 
基因 表达 量 无 变化 ,，58% 基 因 表 达 下 调 ; 而 “TYZ-1 号 "与 正常 浇 水 相 比 , 叶片 中 有 32% 基 因 表达 上 调 , 39% 
基因 表达 量 无 变化 ，29% 基 因 表 达 下 调 ， 根 部 中 36% 基 因 表 达 上 调 ，15% 基 因 表 达 量 无 差异 ，49% 基 因 表 
APU; 在 水 渡 胁 迫 下 ,，“' 佛 奥 ' 根 部 基因 表达 量 与 正常 浇 水 相 比 ，42% 基 因 表 达 上 调 ，9% 基 因 表 达 量 无 变 
化 ,49% 基 因 表 达 下 调 , 叶片 中 30% 的 基因 表达 上 调 , 40% 的 基因 表达 量 不 变 , 30% 的 基因 表达 下 调 ; ‘TYZ-1 
号 ' 在 水 淹 胁 迫 下 与 正常 浇 水 相 比 ，55% 的 基因 在 根部 中 表达 上 调 ，11% 的 基因 在 根部 表达 无 变化 ，44% 的 
基因 在 根部 中 表达 下 调 ， 叶 片 中 29% 的 基因 表达 上 调 ，39% 的 基因 表达 无 变化 ，32% 的 基因 表达 下 调 。 示 
受 胁 迫 时 ，‘TYZ-1 号 ?基因 本 底 表达 普遍 高 于 ' 佛 奥 '( 绿 色 块 多 于 ' 佛 奥 )， 胁 迫 后 "TYZ-1 号 ?基因 表达 量变 
化 不 显著 , 色 块 普遍 较 浅 ， 而 ' 佛 奥 ' 胁 迫 后 深 色 块 多 于 “TYZ-1 号 *， 表 明 ' 佛 奥 基 因 上 调 或 下 调 表达 量变 化 
较 显 著 。 根 据 基因 表达 谱 ， 我 们 从 中 筛选 出 15 个 在 干旱 胁迫 下 表达 量变 化 较 大 的 干旱 胁迫 相关 基因 
(OeAP2-13,0eAP2-33,0eAP2-28,0eAP2-104,0eAP2-82,0eAP2-86,0eAP2-42,0eAP2-23,0eAP2-75,0eAP2-10, 
OeAP2-80,0eAP2-101,0eAP2-79,0eAP2-50,0eAP2-97) 与 15 个 水 淹 胁 迫 下 表达 量变 化 较 大 的 基 
(OeAP2-23,0eAP2-13,0eAP2-101,0eAP2-75,0eAP2-80,0eAP2-12,0eAP2-28,0eAP2-17,0eAP2-50,0eAP2-55, 
OeAP2-104,0eAP2-14,0eAP2-97,0eAP2-2,0eAP2-91) ， 发 现 有 9 个 相同 的 基因 (Oe4P2?-13,Oe4P2-28， 
OeAP2-104,0eAP2-23 ,0eAP2-75,0eAP2-80,0eAP2-101,O0eAP2-50,0eAP 2-97) 4 [rl iil fF BG 7K YAR PI, 这 9 
个 基因 可 能 与 水 分 胁迫 调控 相关 。 
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图 5 油橄榄 AP2/ERF 水 分 胁迫 基因 表达 交互 热 图 
Fig. 5 Interactive heat map of Olea europaea AP2/ERF gene expression under water stress 


2.7 油橄榄 AP2/ERF 基因 qRT-PCR 分 析 
经 特异 引物 tubulinF 和 tubulinR 检测 反 转 录 得 到 油橄榄 cDNA， 采 用 qRT-PCR 分 析 油 橄榄 AP2/ERF 
ft FDL. FDR, ODL. ODR, FSR, FSL, OSR, OSL. FL. FR, OL, OR 中 的 具体 表达 情况 ， 结 果 
示 ( 图 6), OeAP2/ERF 基因 在 FDL, FDR, ODL. ODR, FSR. FSL, OSR, OSL, FL. FR, OL. OR 
表达 量 差异 显著 ， 与 正常 浇 水 处 理 相 比 ，Oe4P2-28 在 ODR 中 表达 量 差异 显著 ，Oe4P2-42 在 FDL PH 
量 差异 显著 ，OeAP2-101、OeAP2-55、OeAP2-42、OeAP2-28 在 FSR 中 均 呈 高 水 平 表达 。 
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不 同 小 写字 母 表示 在 0.05 水 平 上 差异 显著 ， 相 同 小 写字 母 表 示 在 0.05 水 平 上 差异 不 显著 。 
Different lowercase letters mean significant differences at the 0.05 level; Same lowercase letters mean no signific 
ant differences at the 0.05 level. 


图 6 OeAP2/ERF 基因 在 水 胁迫 下 的 表达 情况 
Fig. 6 The expression of OeAP2/ERF gene under water stress 


3 讨论 

用 已 知 功 能 蛋白 作用 机 制 预测 未 知 功能 蛋白 作用 机 制 是 研究 物种 未 知 功能 蛋白 作用 机 制 的 一 种 方法 。 
本 研究 将 油橄榄 AP2/ERF 在 水 胁迫 下 具有 表达 量 的 ERF NE 2S E A 5 38 Xn (Solanum lycopersicum) 
SIERF5(Pan et al., 2012)， 转 基因 大 豆 (Glycine max) GnDREB2A, 7Kf& OsAP37. HRElI. HRE2. JERFI 和 
OsERF109(Oh et al., 2009; Francesco et al., 2010; Zhang et al., 2010; Yu et al., 2017), TU Fir DREB2A(Yoh et 
al., 2006), XVH SpsDRFEB8( 王 雷 等 ，2021) 等 此 类 水 分 调节 功能 基因 共同 构建 蛋白 系统 进化 树 (图 7)， 发 现 
水 胁迫 基因 表达 谱 中 基因 表达 量变 化 较 大 的 5 个 基因 (Oe4P2-75,Oe4P2-97, OeAP2-101,0eAP2-23, 
OeAP2-13) 聚 在 同一 分 支 上 ， 且 该 分 支 中 包含 3 个 已 知 具 有 水 分 调节 功能 的 基因 (HRE1,JERF1,HRE2)， 由 
OeAP2-75,0eAP2-97, OeAP2-101,0eAP2-23, OeAP2-13 Ej HRE1,JERF1,HRE2 具有 相似 调节 功能 ， 

能 与 油橄榄 水 分 调节 密切 相关 。 通过 实时 痰 光 定量 PCR 检测 , 结果 表明 OeAP2-101 基因 在 FSR 及 OSR 
um Le OeAP2-23 在 OSR 中 表达 量 差 异 显著 ，Oe4P2-42 在 FDL 中 表达 量 差 异 显著 ， 水 胁迫 
下 Oe4P2-28 在 2 个 品种 的 叶片 和 根部 表达 量 均 有 显著 性 差异 , Oe4P2-101、OeAP2-55、OeAP2-42、OeAP2-28 
Æ FSR 中 均 呈 高 水 平 表 达 ， 这 进一步 说 明 OeAP2/ERF 参与 油橄榄 水 胁迫 响应 。 
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7 水 分 胁迫 中 具有 表达 量 的 OeAP2/ERF 蛋白 与 其 他 植物 AP2/ERF 水 分 调控 因子 进化 关系 
Fig.7 Evolutionary relationship between OeAP2 /ERF proteins expressed in water stress and AP2/ERF water 


regulatory factors in other plants 


该 研究 中 基因 表达 谱 分 析 发 现 9 个 共同 参与 干旱 与 水 淹 调 节 的 基因 (Oe4P2-73， 
OeAP2-23,0eAP2-75,0eAP2-80,0eAP2-101,0eAP2-50,0eAP2-97), 这 9 FEA CET AK 


OeAP2-28,0eAP2-104, 
榄 水 胁迫 (和 干旱、 水 淹 ) 


下 根 与 叶 的 表达 量 与 非 胁 迫 下 根 与 叶 的 表达 量 有 较 大 的 差异 。 与 非 胁迫 条 件 下 的 基因 表达 量 相 比 , 水 淹 胁 
3À F OeAP2-13 在 ' 佛 奥 ' 中 表达 下 调 ， 在 ‘TYZ-1 号 :中 表达 上 调 ; 水 胁迫 下 的 OeAP2-28 在 ' 佛 奥 ': 和 :TYZ-1 
号 ;的 根 与 叶 中 均 表 达 上 调 ; OeAP2-75 在 干旱 胁迫 下 2 个 品种 的 不 同 组 织 中 均 下 调 表 达 ; OeAP2-101 在 干 


表达 量 都 有 显著 差异 ， 这 些 基 因 在 不 同 组 织 中 表达 量 的 不 同 可 能 就 是 “ 佛 奥 ' 和 “TYZ 


旱 胁迫 下 'TYZ-1 号 ?中 明显 下 调 表达 , 在 水 淹 胁 迫 下 ' 佛 奥 ' 根 部 明显 上 调 表 达 ; OeAP2-104 在 干旱 胁迫 下 2 
个 品种 不 同 组 织 中 均 上 调 ; OeAP2-23, OeAP2-80,0eAP2-50,0eAP2-97 在 水 胁迫 (干旱 、 水 淹 ) 下 不 同 组 织 中 


-1 S pU PESE BS Js 


因 。 转录 组 数据 中 NR 基因 功能 注释 表明 这 9 个 基因 均 为 乙烯 响应 因子 ， 通 过 调控 乙烯 合成 关键 基因 的 表 
达 ， 从 而 调节 乙烯 的 生物 合成 和 信和 号 传导 。Oe4P2?-13,Oe4P2?-75,Oe4P2?-23,Oe4P2?-101,Oe4P2-97 在 油橄榄 
AP2/ERF 蛋白 系统 进化 树 中 属于 ERF B2 亚 族 成 员 ，Maren 等 (2015) 研 究 报道 乙烯 响应 因子 是 激素 和 胁迫 
信号 的 一 个 重要 调节 中 心 , ERF B2 亚 族 成 员 在 低 氧 和 淹 水 反应 中 发 挥 重要 作用 ((Bui et ab, 2015), Xu 等 研 
究 发 现 Sub1A-1(ERF 亚 族 B2 亚 组 成 员 之 一 ) 超 表达 的 同时 促进 乙醇 脱毛 酶 基因 上 调 表 达 ， 抑制 Sub CAKE 


不 耐 受 基因 ) 的 转录 水 平 , 表明 Sub1A-1 在 淹 水 胁迫 应 答 调 控 中 起 关键 作用 (Xu et al., 
橄榄 中 的 表达 模式 与 Sub14-7 相似 ， 在 耐 涝 品种 “TYZ-1 号 "中 表达 上 调 ， 不 耐 涝 品 


2006). OeAP2-13 在 油 
“ 佛 奥 ' 中 表达 下 调 ， 


EW] OeAP2-13 在 “TYZ-1 号 :中 也 通过 促进 乙醇 脱氧 酶 基因 上 调 表达 ,从 而 增强 'TYZ-1 号 耐 涝 性 。 洪林 等 


(2020) 研 究 发 现 RAVs 和 AP2s 调控 植物 的 非 生 物 胁迫 响应 ， 拟 南 芥 RAV1，RAV1L 和 RAV2 在 干旱 胁迫 


下 表达 水 平 下 降 (Fu et al, 2014), Saito 等 (2004) 发 现 过 表达 RAVI 抑制 ABA 降 


fi 3E D (CYP707A1 和 


CYP707A2) 的 活性 ， 从 而 实现 不 依赖 ABA 的 途径 调控 植物 的 非 生物 胁迫 响应 ， 该 研究 中 


OeAP2-28,0eAP2-104, OeAP2-80 在 OeAP2/ERF 蛋白 系统 进化 树 中 属于 RAV, AP2 VI 


族 成 员 ， 在 干旱 胁迫 


下 同样 表达 下 调 ， 推 测 其 通过 促进 或 抑制 干旱 胁迫 相关 酶 的 表达 从 而 调控 油橄榄 干旱 调控 。 

AP2/ERF 转录 因子 在 植物 胁迫 中 调控 机 制 较 复杂 ， 其 启动 子 包含 ABRE( ABA 反应 元 件 结合 蛋白 )、 
DREB( 干 旱 应 答 元 件 结合 蛋白 )、 乙烯 等 顺 式 作 用 元 件 (Zhao etal., 2018; 韩 妙 华 等 , 2020), Sanjana 等 (2016) 
研究 显示 在 植物 胁迫 中 ，DREBs 和 ABREs 等 与 胁迫 相关 的 顺 式 作用 元 件 可 以 和 NAC 转录 因子 结合 从 而 
应 对 植物 胁迫 ，Cao 等 (2020) 研 究 发 现 DREB 结合 WRKY 转录 因子 共同 调节 小 沙 冬青 (4mmopiptanthus 
1a1zzs) 低 温和 王 旱 胁迫 ， 这 一 研究 表明 DREB ABRE 作为 顺 式 作用 元 件 ， 可 同时 调节 AP2/ERF、NAC 和 
WRKY 转录 因子 从 而 应 对 植物 非 生物 胁迫 (Takasaki et al., 2010; WE, 2013). AP2/ERF. WRKY, NAC 
转录 因子 在 植物 水 分 调控 中 扮演 十 分 重要 的 角色 ， 过 表达 TaWRKY 10. VaWRKY14 提高 拟 南 芥 及 转基因 烟 
草 耐 旱 性 (Wang et al., 2017; Zhang et al., 2018); GnDREB2A 参与 大 豆 水 分 亏 缺 调节 (Marinho et al., 2019), 
OsAP37. JERFI. OsERF109 的 过 表达 提高 水 稻 抗 旱 性 (Oh et al., 2009; Zhang et al., 2010; Yu et al., 2017); 
ANACOI9. ANACO55 和 ANACO72 能 使 转基因 植株 抗旱 性 增强 (Tran et al., 2004)， 过 表达 SNACI 增强 转 基 
因 棉 花 (Gossypium spp) 的 耐 旱 性 (Liu et al., 2014); 这 充分 说 明了 AP2/ERF、WRKY、NAC 转录 因子 与 植物 
的 水 胁迫 相关 ， 并 且 在 水 胁迫 中 发 挥 重要 作用 ， 由 此 可 推测 植物 体内 WRKY、NAC 转录 因子 通过 与 
AP2/ERF 中 水 相关 顺 式 作 用 元 件 结合 从 而 共同 调节 水 胁迫 。 课 题 组 前 期 研究 发 现 WRKY 和 NAC 转录 因 
子 参与 油橄榄 水 胁迫 调控 ， 并 且 可 能 与 内 源 激素 ABA 表达 水 平 相关 ， 由 此 可 推测 在 油橄榄 水 胁迫 响应 机 
制 中 ，DREBs、ABREs 作为 顺 式 作用 元 件 ， 同 时 调节 AP2/ERF、WRKY 和 NAC 转录 因子 从 而 共同 调控 
油橄榄 水 胁迫 。 由 于 目前 对 油橄榄 水 胁迫 响应 机 制 及 胁迫 相关 基因 的 综合 研究 还 比较 有 限 , 我 们 对 ‘TYZ-1 
号 "和 人 “ 佛 奥 '2 个 品种 的 油橄榄 进行 了 转录 组 测序 ， 后 续 将 结合 全 基因 组 重 测序 以 揭示 油橄榄 水 胁迫 关键 基 
因 及 潜在 遗传 和 相关 突变 ， 为 油橄榄 分 子 辅助 育种 提供 理论 依据 。 
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